Deeply subrecoil two-dimensional Raman cooling 
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We report the implementation of a two-dimensional Raman cooling scheme using sequential exci- 
tations along the orthogonal axes. Using square pulses, we have cooled a cloud of ultracold Cesium 
atoms down to an RMS velocity spread of 0.39(5) recoil velocities, corresponding to an effective 
temperature of 30 nK (0.15 T rcc ). This technique can be useful to improve cold-atom atomic clocks, 
and is particularly relevant for clocks in microgravity. 
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I. INTRODUCTION 

Types of laser cooling that involve atoms continually 
absorbing and emitting photons cannot in general lead to 
atomic velocity distributions narrower than the recoil ve- 
locity v rec = hk/M where k = 2tt/X is the wavevector of 
photons with wavelength A and M is the atomic mass. By 
contrast, Raman coolingJ3and velocity selective coher- 
ent population trapping @ can reach below this recoil 
limit. These two techniques involve an effective cessation 
of the absorption of the light by the atoms once they have 
reached a sufficiently low velocity. We note that the only 
application of subrecoil cooling of which we are aware Q 
used one-dimensional (ID) Raman cooling. 

Raman cooling has been demonstrated in one, two and 
three dimensions, but deeply subrecoil velocities were 
only obtained in one dimension 0,0. In two and three 
dimensions, the lowest velocity spreads (ID RMS veloci- 
ties) obtained were respectively 0.85 w rcc and 1.34 v rec @- 
Defining the recoil temperature as kBT ICC = Mv^ cc , where 
fee is the Boltzmann constant, these correspond to 0.72 
and 1.80 T lec respectively. In this paper, we report the 
implementation of an efficient 2D Raman cooling scheme 
that has produced velocity spreads as low as 0.39(5) v rcc , 
corresponding to 0.15 T rcc , and that should, under ap- 
propriate circumstances, reach even lower velocities. Our 
technique differs from that used previously in the shape 
of the Raman pulses and the use of sequential excitations 
along the orthogonal axes. 

The use in atomic fountains of ultracold atoms, pro- 
duced by laser cooling in optical molasses, has greatly im- 
proved the accuracy of neutral-atom atomic clocks. Such 
clocks work by launching the atoms vertically through a 
microwave cavity, to which the atoms fall back after a 
Ramsey time as long as about 1 s. The opening in the 
microwave cavity has a diameter typically less than 1 cm, 
so the transverse temperature must be low enough to al- 



low a significant number of the launched atoms to pass 
through the cavity the second time. Atoms that do not 
make it through the second time contribute to the colli- 
sional shift without contributing to the signal. Fountain 
clocks experiments with the coldest atoms achieve RMS 
spreads as low as 2 ?; rec |7j. Under these circumstances, 
many of the atoms are clipped by the second passage 
through the cavity after 1 s of Ramsey time. For signifi- 
cantly longer Ramsey times, for example as envisioned for 
space-borne clocks, even lower temperatures, as obtained 
by subrecoil cooling, will be needed. Note that subrecoil 
longitudinal cooling is not necessarily desirable, because 
the longitudinal thermal expansion of the cloud reduces 
the atomic density and thus reduces the collisional shift. 
For that reason, the present work concentrates on two- 
dimensional Raman cooling of an atomic sample released 
from optical molasses, with a view to providing a valuable 
tool for future atomic clocks. 

The paper is organized as follows: In section |nj we 
summarize Raman cooling theory. In section IIIII we 
present the experimental details, with a stress on the fine 
tuning of the excitation spectrum of the Raman pulses. 
Section IIVI gives the results obtained with our appara- 
tus, in terms of final velocity distribution and cooling 
dynamics. In the conclusion, we summarize and discuss 
our results and their applications. 



II. 



RAMAN COOLING THEORY 
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The theory of one-dimensional Raman cooling in free 
space, described in detail in 0,0, is based on a two-step 
cycle. We consider a cold cloud of cesium atoms initially 
in the hyperfine ground state 6S±/2, F — 3 (at this stage, 
we ignore the Zeeman, mp, degeneracy). First, the atoms 
are placed in the light field of two off-resonant beams 
counter-propagating along the x-axis, with frequencies u>i 
and u>2, and wavevectors ki and k2 such that ki ~ — k2 ~ 
k, where k is a resonant wavevector for the transition 
6S 1/2 -> 6F 3/2 (Fig. [TJ. This light field transfers the 
atoms with some (non zero) velocity along the beams to 
the other hyperfine state 6Si/2,F = 4, while changing 
their velocity by two recoil velocities 7i(ki — k2)/M ~ 
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2Kk/M. The detuning A from the 6P3/2 manifold is 
chosen to be much larger than the hyperfine splitting 
of the upper state and also large enough to avoid one- 
photon excitation. Second, a resonant pulse repumps 
the atoms to F = 3 while giving them the possibility to 
reach an x-velocity close to zero via the emission of a 
spontaneous photon whose cc-component of momentum 
can take any value between ±hk. The Raman detuning 
5, shown in Fig. ^ and defined as S = u>i — u>2 — ^hf, 
is chosen to select atoms with velocities v fulfilling the 
resonance condition 

6 = S LS + 6n+ 45 rcc , (1) 

where 

nf / 1 1 \ nj / 1 1 \ 

LS " 4 V A + ^HF A/ 4 V A -^HF A / ' 

S D = 2k v, 

j- E rcc hk^ 

roc ~ ~hT~2M' 

In these equations, Sli and H2 are the effective electric 
dipole couplings of the Raman beams and whf is the 
hyperfine splitting. The three terms <5ls> <5d and 4£ roc 
are respectively the light shifts, the Raman Doppler ef- 
fect and four times the recoil energy shift. This 3-level 
approach is a good approximation of our problem under 
the following conditions: the polarizations of the Raman 
beams are linear and the detuning A is large compared 
to the upper hyperfine splitting (see the discussion in sec- 
tion Choosing S < selects an initial velocity v 
whose rr-component is opposed to the velocity change, 
which is what we want. Cooling the opposite side of 
the velocity distribution implies repeating the cycle with 




FIG. 1: On top, the Raman transition for the Cesium atom. 
At the bottom, the beam configuration for a velocity selective 
transition. We have shown the usual conditions where 5 and 
A are negative. 



the directions of ki and k2 reversed from that shown in 
Fig.ffl 

Because of common mode rejection, only relative fre- 
quency noise between the Raman laser beams affects the 
Raman selectivity. By phase locking them relative to 
each other, this difference- frequency noise can be made 
much smaller than the noise of their separate frequencies, 
and negligible. The excitation spectrum is then fully de- 
termined by the shape and amplitude of the pulses. A 
careful tailoring of the pulse shape allows a precise exci- 
tation spectrum that does not excite atoms with a zero 
velocity along the x-axis (see for example Fig. 0^,). By 
repeating the cooling cycle a large number of times, one 
forces the atoms to perform a random walk in velocity 
space until they hit the zero velocity state, a so-called 
dark state, where they tend to accumulate. 

Our two-dimensional Raman cooling is a direct exten- 
sion of the one-dimensional case, where the cooling cycles 
are alternatively applied to the x and the y directions. 

The first Raman cooling experiments used Black- 
man pulses, which feature a power spectrum with very 
small wings outside the central peak, hence reducing off- 
resonant excitations. Although this might seem to be a 
very desirable feature, later work |5| showed experimen- 
tally and theoretically that square pulses, which produce 
an excitation spectrum featuring significant side lobes 
and a discrete set of zeros, give a better cooling in the 
one-dimensional case. It is also expected to be better in 
the two-dimensional case 0, Q . The dynamics of Raman 
cooling, as well as that of VSCPT, are related to non 
Gaussian statistics called Levy flights 0, 0] . More pre- 
cisely, for an excitation spectrum varying as v a around 
v = 0, the width of the velocity distribution scales with 
the cooling time 8 as _1 / Q . However the atoms effi- 
ciently accumulate in the cold peak of the distribution 
when — * 00 only if a is greater than or equal to the 
dimensionality of the problem. Square pulses, for which 
a = 2, appear to be suitable for 2D cooling, and the 
present work concentrates on them. 



III. EXPERIMENTAL SETUP 

A. Laser system 

Raman cooling requires two laser beams whose fre- 
quency difference is locked to a frequency close to the 
hyperfine frequency whf/27t ~ 9.2 GHz. The most com- 
mon methods used to generate the two frequencies in- 
clude direct electronic phase locking of two free running 
lasers 0, acousto-optic modulation 0] > an d electro- 
optic modulation p| . We used a different a ppr oach based 
on current modulation of a laser diode |l3l Il4| , as shown 
in Fig. [21 An extended-cavity master diode laser at 
852 nm, with a free spectral range of 4.6 GHz, is current- 
modulated at o>hf /2 in order to generate sidebands sep- 
arated by whf ■ The fraction of the power in the two first 
order sidebands, measured with an optical spectrum an- 
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alyzer, is about 50%. The carrier is filtered out with a 
solid etalon having a free spectral range of 9.2 GHz and a 
finesse of 8, and the remaining beam is used to injection- 
lock two slave diodes. The slave currents are adjusted in 
order to lock one slave to one sideband and the other slave 
to the other sideband. In the spectra of the slaves, the 
total contamination from the carrier and any of the un- 
wanted sidebands is less than 1% of the total power. The 
phase coherence of the sidebands is fully transfered onto 
the slaves and the beatnote spectrum of the two slaves is 
measured |25j to be 1 Hz wide. This includes contribu- 
tions from the linewidth of the microwave generator used 
to modulate the master laser, the mechanical vibrations 
of the laser and optical system (but not of those mirrors 
after the fibers), and the resolution of the measurement 
apparatus. After transport in optical fibers, 40 mW are 
available in each Raman beam. 

As pointed out previously, cooling of opposite sides of 
the velocity distribution requires interchanging the di- 
rections of ki and k2. This is done by interchanging 
the injection currents of the slaves so that the sidebands 
to which they lock are interchanged, thus interchanging 
their roles |l5j. The switching time, measured by mon- 
itoring the transmission of each Raman beam through 



,4.6 GHz 




FIG. 2: Schematic of the setup used to generate the Raman 
beams, showing the optical spectrum of the beams at vari- 
ous places. AOM: acousto-optic modulator; PBS: polarizing 
beam-splitter; PC: Pockels cell. The AOMs before the slave 
lasers are frequency shifters, and the AOM between the PBSs 
is a shutter. 



a confocal cavity [l^, is found to be about 30 fis for a 
complete switch, similar to what was observed in Ref. Il5l 
In the cooling experiment, we allow an extra 20 fj,s as a 
safety margin. The swapping of the beams between the x 
and the y directions is done with Pockels cells and polar- 
izing beam splitters (Fig. in less than a microsecond. 
The extinction ratio of the Pockels cell switches is about 
100. 



B. Experimental details 

A magneto-optical trap inside a glass cell is loaded 
with a few times 10 7 atoms from a chirped-slowed atomic 
beam. The cloud has an RMS width of about 1 mm. Af- 
ter additional, 70 ms-long, molasses cooling, the atoms 
are dropped, pumped into F = 3, and Raman cooled 
for 25 ms, before they fall out of the Raman beams. As 
shown in Fig. [21 the Raman beams are in the horizontal 
plane, along the x and y axes, providing cooling perpen- 
dicularly to the vertical direction. We found that con- 
trolling the horizontality of the Raman beams at a level 
of a few thousandths of a radian is enough to ensure 
that gravity does not perturb the cooling, but an error 
as large as 0.01 rad has a noticeable effect. The waist 
of the beams (radius at 1/e 2 of peak intensity) is 4 mm 
and they all have nominally the same power. When the 
atoms are dropped, they are slightly above the center of 
the beams, and after 25 ms of cooling, they are at an 
approximately symmetric position below. 

The repumping is provided by a retro-reflected verti- 
cal beam, tuned to the F — 4 — > F 1 = 3 transition, 
with an intensity a few times the saturation intensity and 
with the reflected polarization rotated in order to avoid 
a standing wave effect. The momentum of the photons 
absorbed from the repumping beams has no effect on the 
transverse velocity, and leads only to momentum diffu- 
sion in the vertical direction. 

The velocity distribution along x or y is measured by 
Raman spectroscopy |l6j . i.e. by transferring a narrow 
velocity class from F — 3 to F — 4 with a long Ra- 
man 7r pulse. Two centimeters below the Raman beams, 
the atoms fall through a sheet of light tuned to the 
F = 4 — > F 1 = 5 cycling transition. The integrated 
fluorescence collected by a photodiode is proportional to 
the number of atoms transfered to the F — 4 state. By 
scanning the Raman detuning S for a succession of iden- 
tically cooled atomic samples, one can probe all the ve- 
locities and reconstruct the velocity distribution. There 
is a small, uniform background signal; after subtraction 
of this background, we obtain velocity distributions such 
as that shown in Fig. [5] 

Raman cooling as described here is essentially a 3-level 
scheme and our experiment requires the Zeeman sub- 
levels to be degenerate within each hyperfine level. Good 
subrecoil Raman cooling can be achieved only if any Zee- 
man splitting is small compared to the Raman Doppler 
shift associated with a single recoil velocity, which is 
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Raman beams 




sheet of resonant light 



FIG. 3: Geometry of the experiment. The Raman beams 
are in the horizontal plane, while the repumper beams are 
vertical. The polarizations are linear along the directions in- 
dicated by the arrows. The MOT beams are not shown. The 
atoms in F = 4 are detected by fluorescence when they fall 
through a horizontal sheet of light resonant with the transi- 
tion F — 4 — > F' — 5, located 2 cm under the MOT position. 



8.2 kHz. This is ensured by reducing the DC stray mag- 
netic field with an opposing applied external magnetic 
field, and further reducing the DC and AC residual fields 
with a /x-metal shield. Raman spectroscopy with non- 
velocity-selective, co-propagating beams and long pulses 
(300 /is) is used to optimize the field zeroing by adjusting 
for minimum spectral width. The Raman spectrum has 
a full width half maximum (FWHM) of 0.5 kHz, corre- 
sponding to a residual stray field smaller than 100 /iG, 
and equivalent to the Raman Doppler shift of atoms with 
a velocity w r cc/16. 



C. Excitation spectrum 

The polarizations of each pair of Raman beams are 
crossed-linear in order to ensure that, because the de- 
tuning |A| is large compared to the hyperfine splitting 
(600 MHz) of the excited state, the light shifts are nearly 
the same for all the Zeeman sublevels of the ground 
state [lTj • Under those conditions, the effective electri- 
cal dipole coupling SI corresponding to an intensity / has 
a value O = r^/0.67 I/2I , where T = 2w X 5.2 MHz 
is the natural lincwidth of the excited state, and Iq = 
1.1 mW/cm 2 is the saturation intensity for the strongest 
transition. The Raman detuning 6 has to be negative to 



cool the atoms, and is chosen in such a way that atoms 
with a zero velocity are resonant with the first zero point 
of the excitation spectrum [2f|. We extend Eq. Q by 
defining the effective detuning seen by these zero veloc- 
ity atoms as 



5 a = 5 - 6_ 



LS 



4<5 r , 



so that 



r 2 / 



<5o = 5 + 0.67 

8i V A - W HF 



1 



A + w H f 



2hk 2 
M 



(2) 



The excitation spectrum, defined as the probability V of 
undergoing a Raman transition for any atom seeing an 
effective Raman detuning <5 CX , is given for a square pulse 
of length t by the Rabi formula: 



n 2 



S 2 



I J dr. 



(3) 



where SIr ~ ■ is the Raman Rabi frequency. The 

coefficient C depends on the initial state in the F = 3 
manifold. It has a mean value C = 0.28 and a total 
spread of ±20%. The value of 6 defined by Eq. © 
must fulfill V(5q) = 0. That is, for zero velocity atoms, 
t^JS 2 + Sip = 2-7T. In the above, we have assumed that 
S7i — f^2 = SI- In fact, the intensities 1\ and I2 of the 
Raman beams may differ by as much as 20%. To take 
this into account, one would have to write Eq. JSJ) in the 
form of Eq. (QJ , and replace / in the definition of S1r with 
sfhh- 

The detuning 5 depends on the light intensity in two 
different ways. Firstly, there is the differential light shift 
#ls between the two hyperfine levels, which is propor- 
tional to the light intensity /. Secondly, according to 
Eq. (PJ, the frequency of the first zero of the excita- 
tion spectrum changes due to a "saturation" effect, as 
soon as S7r£ is not small compared to 1, that is to say 
when the maximum transfer probability is not small com- 
pared to unity. These dependences on the light intensity 
lead to complications because the Raman beams have 
a Gaussian profile and are never perfectly spatially ho- 
mogeneous. Depending on their position in the beams, 
different atoms have different resonance conditions. In 
previous experiments 0, 0], the differential light shift 
was reduced by using a detuning A large compare to the 
hyperfine splitting cjhf , and the saturation effect was re- 
duced by setting the maximum transfer probability to 
0.5, thus using ir/2 pulses instead of ir pulses. 

Because of a limitation in the available laser power, we 
worked at a detuning | A| of only 20 GHz, with pulses hav- 
ing a maximum transfer efficiency of 80% (0.77r-pulses). 
As shown below, although not negligible, the saturation 
shift and the light shift partially cancel out, thus limiting 
unwanted excitation of the dark state. 

We measure the excitation spectrum by determining 
the transfer efficiency of such Raman pulses as a function 
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of the Raman detuning 6, with non-velocity-selective co- 
propagating beams. Figure 0] shows the excitation spec- 
trum when the atomic cloud is centered on the Raman 
beams (a), resulting in a fairly homogeneous illumina- 
tion, and when the cloud is on the edge of the beams 
(b), resulting in an inhomogeneous illumination. The 
arrows show the position of the Raman detuning corre- 
sponding to the dark state. Because the light shift and 
the saturation respectively shift and broaden the spec- 
trum, locations in the cloud exposed to different light 
intensities yield different positions and shapes of the ex- 
citation spectrum. The spectra of the individual atoms 
contribute inhomogeneously to the measured spectrum. 
In both case (a) and (b), it appears that the spectrum 
is not fully symmetric and that the first minimum on 
the positive detuning side does not go as close to zero as 
does the minimum corresponding to the dark state. This 
comes from the fact that the saturation and the light 
shift have opposite effects on the position of the zero 
on one side of the spectrum, and similar effects on the 
other side. This effect becomes very visible in the inho- 
mogeneous illumination case. However, the cancellation 




Effective Raman detuning § ex (kHz) 

FIG. 4: Excitation spectrum of 50 /is-long Raman pulses, 
with a maximum transfer probability of 0.8. In figures (a) 
and (c), the atomic cloud is centered on the Raman beams 
(nearly homogeneous illumination), in figures (b) and (d), 
the cloud is on the side of the beams (inhomogeneous illu- 
mination). The position of the dark state is indicated by the 
arrows. The dashed line in figure (a) is the parabola that 
matches the second derivative of the spectrum at the position 
of the dark state. 



is not perfect, and the dispersion of the Raman coupling 
through the dependence of C on the Zeeman sublevel also 
leads to a "blurring" of the spectrum. As a result, the 
dark state features a small excitation probability, even in 
the homogeneous illumination case. On the time scale of 
our cooling sequence, this has not proven to be of impor- 
tance. 

With a detuning |A| = 20 GHz, our 50 /is-long Raman 
0.7-7T pulse corresponds at the center of the beam to a 
mean Raman Rabi frequency £!r = Cr 2 /4|A| ■ I/Iq = 
2ir x 7 kHz, and an intensity / ~ 80 mW/cm 2 for each 
beam. The probability of one-photon excitation is of the 
order of r/| A| per pulse. It results in a total probability 
of excitation of 10% for a typical cooling sequence of 300 
pulses, which is low enough to avoid any significant per- 
turbation of the cooling. The choice of a pulse length of 
50 /is means that the excitation spectrum covers most of 
the initial velocity distribution, as seen on Fig. [SJ 



IV. RESULTS 

The elementary Raman cooling cycle is a 50 /is square 
Raman pulse with the lo\ beam along some direction, e.g. 
+x (as in Fig. ^) , and the u>2 beam along the opposite di- 
rection, e.g. —x, followed by a 10 fis resonant repumping 
pulse from counter-propagating beams along z. In this 
example, the elementary cycle provides cooling along the 
+x axis. The rise and fall times of the Raman pulse are 
less than 1 fj,s. In our experiment, the elementary cool- 
ing cycles are applied in pairs along a given direction. 
A complete cycle consists of 4 pairs of elementary cool- 
ing cycles applied successively along the directions +x, 
+y, —y and — x. The full cooling sequence is typically 
composed of 40 complete cycles. 

As pointed out previously, the switching time between 
perpendicular directions, limited by Pockels cell switch- 
ing, is instantaneous with respect to the experimental 
timescale, but the switching time between parallel direc- 
tions is 50 /is. Using pairs of identical elementary cycles 
reduces the total number of complete cycles, and thus 
reduces the total time spent switching the beams. There 
is almost no loss of efficiency resulting from the applica- 
tion of two successive elementary cycles along the same 
direction because after only a few complete cooling cy- 
cles, the velocity distribution becomes narrow enough so 
that the atoms are in the tail of the excitation spectrum 
and the excitation probability per Raman pulse is small 
compared to 1 for most of the atoms. 

Immediately after the cooling, the velocity distribu- 
tion is probed with a 500 fis square ir pulse applied 
along the x or the y axis. Figure |SJd shows the veloc- 
ity distribution measured along the x axis after 320 ele- 
mentary cycles (40 complete cycles), corresponding to a 
22.4 ms cooling seque nce. The velocity spread, defined as 
a = FWHM/ \/8 In 2, is reduced from about 4 v rcc after 
molasses cooling, corresponding to an effective temper- 
ature of 3 fiK, to 0.39(5) u ra c |27j . corresponding to an 
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effective temperature of 30 nK. Note that a 20% con- 
tribution from the probe excitation linewidth is removed 
by deconvolution. A similar velocity distribution is mea- 
sured along the y axis. For either of these ID velocity 
distributions, about 80% of the area is contained in the 
central cold peak. Most of the rest of the atoms are con- 
tained in two secondary peaks at positions where both 
excitation spectra are close to zero, as seen in Fig. [SJi. 
From this we conclude that in 2 dimensions, about 
60% of the atoms are in the 2D cold peak. 

It is worth noting that because we operate in the case 
where a is equal to the dimensionality of the cooling, 
the number of atoms in the cold peak of the 2D velocity 
distribution is predicted 0, |t| not to change any more as 
soon as we enter the subrecoil regime and the cold peak 
is clearly separated from the rest of the distribution. The 
main effect of continued cooling is to reduce the width 
of the cold peak without changing the number of atoms 
it contains. This is indeed what we observe as shown in 
Fig. El 

During the entire Raman cooling process, the 2D veloc- 
ity width is reduced by a factor of 10, while the size of the 
cloud along x and y does not change significantly, leading 
to a 100-fold increase of the 2D phase space density. Of 
course, taking into account the heating and expansion 
along the z direction (which we did not measure), the 
increase in the 3D phase space density would be smaller, 
perhaps a few tens. We note however that Raman cooling 
in free space is not a promising way to increase the phase 
space density to reach quantum degeneracy because the 
spatial density stays roughly constant. This means that a 
great deal of additional cooling is required. On the other 
hand, a high spatial density is not desirable for clock ap- 
plications, because it leads to collisional frequency shifts. 

An important issue is the isotropy of the velocity distri- 
bution in the cooling plane. The cooling scheme is funda- 
mentally anisotropic, and we only measured the velocity 
spread along the cooling directions x and y. However 
there are good reasons to believe that the peak of cold 
atoms at the center of the final distribution is isotropic 
(in two dimensions). As indicated in [T(ij . the dynamics 
of each individual atom is dominated by fairly distinct 
phases, where it either performs a random walk in ve- 
locity space outside the subrecoil range, or stays close to 
the dark state, in the subrecoil range, until it gets ex- 
cited and resumes the random walk. The argument for a 
final isotropic velocity distribution relies on two consid- 
erations. Firstly, the total excitation probability for an 
atom close to the dark state during a complete cycle is 
isotropic. Indeed, close to the null velocity, the excitation 
probabilities for an elementary cycle along ±x or ±y are 
small compared to 1 and proportional respectively to v\ 
and Vy (power law with a coefficient a = 2), leading to 
a total excitation probability for the complete cycle pro- 
portional to the sum v% = v% + v 2 y . This probability only 
depends on the "distance" v r from the dark state, and 
is therefore isotropic. Secondly, atoms excited from the 
subrecoil range perform several steps during their ran- 



dom walk before going back to the subrecoil range at 
a random point which is uncorrelated with the previous 
position they occupied close to the dark state. 

The combination of the effectively isotropic excitation 
of the subrecoil atoms and the homogeneous filling of the 
subrecoil region should produce an isotropic cold peak. 
We checked that a simple Monte Carlo simulation ignor- 
ing Zeeman sublevels gives a perfectly isotropic distribu- 
tion in the x — y plane. 

We also studied the experimentally measured velocity 
spread a as a function of the cooling time 9. The results 
are shown on Fig. As stated in section II I II subre- 
coil cooling theory predicts that the velocity spread is 
described at long times by the power law -1 / Q , where 
a is the excitation spectrum power law coefficient. Our 
data do not cover a range of cooling times large enough 
to fully enter the asymptotic regime, and to allow an ac- 
curate experimental determination of a. However, Fig.|S| 




Atomic velocity (v. ) 




FIG. 5: (a) Velocity distribution along the x axis before Ra- 
man cooling (circles), and after 15 ms of 2D Raman cooling 
(heavy line), (b) Velocity distribution along the x axis after a 
22.4 ms 2D Raman cooling sequence consisting of 320 Raman 
pulses. Also shown in both graphs in dotted and dashed lines 
are the excitation spectra used to cool along the directions x 
and —x, tuned to match their first zeros on the dark state. 
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Cooling time (ms) 



FIG. 6: Filled squares: Velocity spread a as a function of the 
cooling time O. The straight line is the power law /30 _1//2 , 
whose multiplicative factor j3 has been adjusted manually to 
make the line fit the data at long times. Open circles: Fraction 
of atoms in the cold peak determined as explained in the text. 



shows that our data is compatible with cooling dynamics 
described by the theoretical power law O -1 / 2 (a = 2) at 
times longer than 10 ms. 

It is expected that because of the experimental im- 
perfections, the velocity spread would eventually reach 
a finite value at long times. Nonetheless, Fig. does 
not show any evident saturation, which indicates that a 
longer cooling time, for instance in micro-gravity, would 
lead to an even smaller velocity spread. 

One of the potential limitations of Raman cooling is 
photon reabsorption fl9| . However, with a spatial density 
less than 10 10 cm -3 and an optical depth of the order of 
unity, we do not expect photon reabsorption to hinder 
the cooling [20|. In clock applications, the density would 
typically be significantly lower than in these experiments 
and reabsorption should be completely negligible. 

To keep the cooling sequence simple, we use Raman 
pulses with a fixed length. The choice of 50 /is-long pulses 
is convenient because the resulting excitation spectrum 
covers most of the initial velocity distribution, and more 
importantly, covers the maximum excursion range of the 
atoms during their random walk. Indeed, an atom close 
to zero velocity can be pushed away from the center of 
the velocity distribution by a maximum amount of about 
4 w rcc : two recoil velocities during the Raman transition 
plus one or two (or exceptionally more) during the re- 
pumping process. The excitation spectrum covers 5 v Iec 
between the two first zeros. 

It should be possible to use longer pulses (narrower 
excitation spectrum) in combination with short pulses 
(wider excitation spectrum) that recycle atoms far from 
zero velocity The increased filtering effect of longer 
pulses produces an narrower distribution at the cost of a 
smaller number of atoms in the cold peak. In any case, 
the best cooling strategy results from a tradeoff between 
the width of the distribution and the fraction of atoms 



in the cold peak, and depends on the total cooling time 
available. 



V. CONCLUSIONS 

Our scheme produces a narrower velocity distribution 
with respect to the recoil velocity than what was previ- 
ously achieved with 2D Raman cooling [||. There are 
two main differences from that previous scheme, where 
the cooling was performed in a vertical plane, from the 
four directions at the same time. First, we use one pair of 
Raman beams at a time, in order to avoid unwanted ex- 
citation of the dark state. Indeed, in Ref. |fg, the careful 
use of circular polarization avoided diffracting the atoms 
from standing waves created by counter-propagating 
beams of the same frequency, but higher order pho- 
ton transitions of the type (wi , k) (w 2 , — k) (u)2 , k) (uji , — k) 
would still be able to transfer 4 v rcc to the atoms. Second, 
since we cool in the horizontal plane, gravity has no effect 
on the velocity components which are being cooled and 
does not accelerate atoms out of the cold peak, allowing 
for a more effective cooling. 

To be used in an cold-atom atomic clock, Raman cool- 
ing has to be coupled with a launching mechanism like 
moving molasses. An easy solution is to first launch the 
atoms and then collimate them with Raman cooling on 
their way towards the first microwave cavity. The cooling 
time depends on the size of the Raman beams and the 
launch speed. 

In that perspective, our setup performs quite well in 
comparison with a newer scheme relying on sideband 
cooling in optical lattices [2l]], which has produced po- 
larized samples with a velocity spread of 0.85 u rec in 3D, 
in a fountain- like geometry |22( . While it is appealing for 
its simplicity, sideband cooling has a fundamental limit 
for the lowest velocity spread achievable, which is about 
0.7 times the recoil velocity associated with the wavevec- 
tor of the lattice used to trap the atoms [23j. Raman 
cooling has no such a limitation. 

It is worth noticing that, although deeply subrecoil ve- 
locities are obtained for long cooling times, only accessi- 
ble in micro-gravity, 2D Raman cooling can still provide 
subrecoil velocities in a few milliseconds, as shown in 
Fig. |SJ Implementing the scheme on a moving-molasses 
earth-bound fountain, where a 1.5 cm- interaction region 
with the Raman beams combined with a typical launch 
velocity of 5 m/s leads to an interaction time of 3 ms, 
would yield a substantial improvement in terms of bright- 
ness of the atomic source, reducing the transverse veloc- 
ity spread from a few recoils to less than a recoil velocity. 

For a micro-gravity-operated atomic beam, the im- 
provement would be even more dramatic because the 
launch velocity can be much smaller than in a fountain, 
making the interaction time with the Raman beams much 
longer. The maximum cooling time is more likely to be 
limited by the maximum longitudinal heating acceptable. 
How the increase of brightness translates into an increase 
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of the stability of a space-borne atomic clock depends 
on geometrical details and on the factors which actually 
limit the stability and/or the accuracy. For simplicity, 
let us assume that the atomic cloud is severely clipped 
by the opening of the second microwave cavity, as is the 
case in current fountain clocks, and that the signal-to- 
noise ratio is the main limiting factor of the short-term 
stability (2^. Reducing the transverse velocity spread 
from 2 w r oc to 0.4 v xec (our current result) would increase 
the flux of atoms through the cavity by a factor of 25. 
That would translate into a 5-fold increase of the signal- 
to-noise ratio, leading to a 5-fold increase of the short- 
term stability for a given averaging time, or a 25-fold re- 



duction in the averaging time needed to achieve a given 
short-term stability. 



Acknowledgments 

We thank F. Bardou and C. Ekstrom for very helpful 
discussions. We also thank C. Ekstrom, W.M. Golding 
and S. Ghezali for early contributions to the experimental 
apparatus. This work was funded in part by ONR and 
NASA. 



[1] M. Kasevich and S. Chu, Phys. Rev. Lett. 69, 1741 
(1992). 

[2] A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, 
and C. Cohen- Tannoudji, Phys. Rev. Lett. 61, 826 

(1988) . 

[3] A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, 
and C. Cohen- Tannoudji, J. Opt. Soc. Am. 6, 2112 

(1989) . 

[4] M. BenDahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996). 
[5] J. Reichel et al, Phys. Rev. Lett. 75, 4575 (1995). 
[6] N. Davidson, H.-J. Lee, M. Kasevich, and S. Chu, Phys. 
Rev. Lett. 72, 3158 (1994), note that this reference 
quotes 2D and 3D RMS velocities which we converted 
to ID RMS velocities by dividing by \fl and \/3 respec- 
tively. Note also that our definition of T rcc is different 
from that of Davidson et al. 
[7] S. Jefferts, T. Heavner, E. Donley, J. Shirley, and 
T. Parker, in Proc. 2003 Joint Mtg. IEEE Intl. Freq. 
Cont. Symp. and EFTF Conf. (2003), p. 1084. 
[8] K. Moler, D. S. Weiss, M. Kasevich, and S. Chu, Phys. 

Rev. A 45, 342 (1992). 
[9] F. Bardou, J.-P. Bouchaud, A. Aspect, and C. Cohen- 
Tannoudji, Levy statistics and laser cooling (Cambridge 
University Press, 2001). 

[10] F. Bardou, J. P. Bouchaud, O. Emile, A. Aspect, and 
C. Cohen-Tannoudji, Phys. Rev. Lett. 72, 203 (1994). 

[11] G. Santarelli, A. Clairon, S. N. Lea, and G. Tino, Opt. 
Commun. 104, 339 (1994). 

[12] P. Bouyer, T. L. Gustavson, K. G. Haritos, and M. A. 
Kasevich, Opt. Lett. 18, 649 (1993). 

[13] J. Ringot, Y. Lecoq, J. C. Garreau, and P. Szriftgiser, 
Eur. Phys. J. D 7, 285 (1999). 

[14] L. Goldberg, H. F. Taylor, J. F. Weller, and D. Boom, 
Electron. Lett. 19, 491 (1983). 

[15] K. Szymaniec, S. Ghezali, L. Cognet, and A. Clairon, 
Opt. Commun. 144, 50 (1997). 

[16] M. Kasevich, D. S. Weiss, E. Riis, K. Moler, S. Kasapi, 
and S. Chu, Phys. Rev. Lett. 66, 2297 (1991). 

[17] J. D. Miller, R. A. Cline, and D. J. Heinzen, Phys. Rev. 
A 47, R4567 (1993). 

[18] J. Reichel, Ph.D. thesis, University of Paris VI (1996). 

[19] Y. Castin, J. I. Cirac, and M. Lewenstein, Phys. Rev. 
Lett. 80, 5305 (1998). 

[20] H. Perrin, A. Kuhn, I. Bouchoule, T. Pfau, and C. Sa- 



lomon, Europhys. Lett. 46, 141 (1999). 
[21] S. E. Hamann, D. L. Haycock, G. Klose, P. H. Pax, I. H. 
Deutsch, and P. S. Jessen, Phys. Rev. Lett. 80, 4149 
(1998). 

[22] P. Treutlein, K. Y. Chung, and S. Chu, Phys. Rev. A 63, 
051401 (2001). 

[23] A. Kastberg, W. D. Phillips, S. L. Rolston, R. J. C. 
Spreeuw, and P. S. Jessen, Phys. Rev. Lett. 74, 1542 
(1995). 

[24] F. P. DosSantos, H. Marion, S. Bize, Y. Sortais, A. Cla- 
iron, and C. Salomon, Phys. Rev. Lett. 89, 233004 
(2002). 

[25] The linewidth of the beatnote spectrum was measured 
by recording the beating of the two beams on a fast pho- 
todiode. The photo-signal was mixed with the signal of 
an auxiliary microwave generator tuned to a frequency 
close to 9.2 GHz. The mixing signal was analyzed with 
an FFT spectrum analyzer. 

[26] The Raman detuning S is experimentally adjusted (by 
optimizing the final velocity distribution) in such a way 
that the velocity class resonant with the first zero point 
of the excitation spectrum is the same when we cool both 
sides of the velocity distribution. However, because we do 
not know precisely the value of the light shifts, nor do we 
have a perfect calibration of S, such a dark state is not 
necessarily the zero velocity state. In fact, this degree of 
freedom can be used to tune the direction of propagation 
of the atoms after the cooling. 

[27] Uncertainties quoted in this paper represent one standard 
deviation, combined statistical and systematic uncertain- 
ties. 

[28] Simulations of the cooling process show that the 2D ve- 
locity distribution consists of one peak centered on the 
zero velocity (v x — 0, v y =0), and 4 additional cold 
peaks centered approximately at positions (v x = ±5u re c, 
v y = 0) and (v x — 0, v y — ±5u rC c). These 4 peaks cor- 
respond to secondary dark states and are experimentally 
determined to contain about 10% of the atoms each. 

[29] One could choose instead to take advantage of the im- 
proved collimation by reducing the initial number of 
atoms launched in order to reduce the collisional shift, 
which is a major source of inaccuracy in laser cooled Ce- 
sium atomic clocks. However, recent developments |24[ 
suggest that the collisional shift can be accurately mea- 
sured and accounted for. 



